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CL|' Abstract 

We investigate the power-suppressed corrections to structure functions in flavour 
singlet deep inelastic lepton scattering, to complement the previous results for 
^ I the non-singlet contribution. Our method is a dispersive approach based on an 

' analysis of Feynman graphs containing massive gluons; and our results agree 

with those obtained from leading infrared renormalon contributions. As in 
non-singlet deep inelastic scattering we find that the leading corrections are 
proportional to We find that the singlet contribution becomes important 

at small x. 
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1 Introduction 



The study of structure functions in deep inelastic lepton scattering (DIS) has received a 
great impetus from the increasing quantity and kinematic range of the HERA data |l|]. 
Although these functions cannot be calculated using perturbative QCD, their asymptotic 
scaling violations (logarithmic dependence) can be predicted and used to measure the 
strong coupling a^. (The same is true of the study of final-state properties, such as frag- 
mentation functions and event shape variables.) 

One problem with the measurement of as using scaling violation, in either structure or 
fragmentation functions, is that there is dependence associated with power-suppressed 
(higher- twist) contributions, in addition to the dominant logarithmic dependence. These 
contributions need to be estimated in order to make use of the wide coverage of HERA. 

Recently so-called 'renormalon' or 'dispersive' methods of estimating power-suppressed 
terms have been applied to a wide variety of observables. By looking at the behaviour 
of the QCD perturbation series in high orders, one can identify unsummable, factorially 
divergent sets of contributions (infrared renormalons |§, i|, H]) which indicate that non- 
perturbative power-suppressed corrections must be included. The Q^-dependence of the 
leading correction to a given quantity can be inferred, and by making further universality 
assumptions one may also estimate its magnitude. Tests of these ideas provide information 
on the transition from the perturbative to the non-perturbative regime in QCD. In partic- 
ular, one can investigate the possibility that an approximately universal low-energy form 
of the strong coupling may be a useful phenomenological concept p, |^ . 

Such an approach has been applied with some success to e+e~ fragmentation functions 
and event shape variables 0, 10, 11], and to flavour non-singlet DIS structure functions 



1^, 0, ITSf, fragmentation functions [T^ and event shape variables Note that special 



care has to be taken when considering non-inclusive observables such as fragmentation 



functions and event shape variables |16]. There are also renormalon model results for 



photon-photon scattering [|I^ and structure functions in flavour singlet DIS |T8|. 

In the present paper we perform the singlet structure function calculation using the 
dispersive method. We flnd agreement with the renormalon model result: as in non-singlet 
DIS, the predicted leading power corrections to these quantities are proportional to 1/Q^, 
but their functional forms are different. In particular there are contributions proportional 
to (logQ2)/Q2^ which are not found in the non-singlet case. The hypothesis that they are 
related to a universal low-energy strong coupling implies that their magnitudes are given 
by a universal non-perturbative parameters. 

In the following section we review the approach of reference . Section § presents the 
standard leading-order perturbative treatment of DIS structure functions. In section ^ we 
estimate the power-suppressed corrections using the method outlined in section |[ Our 
results are summarized briefly in section ^. 

2 The Dispersive Approach to Power Corrections 

We assume that the QCD running coupling a;s(fc^) can be deflned for all positive fc^, and 
that apart from a branch cut along the negative real axis there are no singularities in the 



1 



complex plane. It follows that we may write the dispersion relation: 

«.(^^)=- r^xT^Ps(/i^), (2.1) 

where the 'spectral function' ps represents the discontinuity across the cut: 

Ps(/i') = T^lasip'en - asifi^e-'n] = T^Disc a,( V) • (2.2) 

We now consider the calculation of some observable F in an improved one-loop approxi- 
mation which takes into account one-gluon contributions plus those higher-order terms that 
lead to the running of a^. As discussed in |Q, we expect, for squared Feynman diagrams 
containing a single gluon, that 

where the characteristic function J-'{fi'^ /Q'^) is obtained by one-loop evaluation of F (divided 
by as) with the gluon mass set equal to /i 0, ||. The first term on the right-hand side 
represents the contributions in which a single gluon is produced or exchanged, while the 
second represents those with more complex final or virtual states (e.g. the 'decay products' 
of a virtual gluon, which contribute to the running of a^). 

We can eliminate 0:^(0) by means of the dispersion relation ( p.l|) : 



F = as{0)m + r ^ Psif^') ^(pVQ') , (2.3) 
Jo 



Jo ^ J 



By rotating the integration contour separately in the two terms of the discontinuity, we 
obtain the following contribution to the observable F: 

F = as{p') g{pVQ^) + /ooN , (2.5) 

Jo 



where, setting p? IQ"^ = e. 



g{e) = --^Disc ^(-e) , (2.6) 

ZTTZ 

and the functional Joo is some integral around a contour at infinity. 

Non-perturbative effects at long distances are expected to give rise to a modification 
in the strong coupling at low scales, Sasip"^) = Ois{p'^) — ^^"'"(/i^), a^^i^p^) being the 
perturbatively-calculated running coupling. This generates the following non-perturbative 
correction to the perturbative prediction: 

SF = 6as{p') QipVQ') . (2.7) 

Jo p^ 

Since Sas{p'^) is limited to low values of p"^, there is no contribution to SF from the 
contour at infinity. Furthermore, the asymptotic behaviour of 6F at large is controlled 
by the behaviour of ^(e) as e — 0. We see from ( |2.(j| ) that no terms analytic at e = 
can contribute to SF. On the other hand non-analytic terms at small e do contribute, the 
relevant terms being: 

:F^a,^V-e =^ 5F = -^^, (2.8) 

ZTT 71 Q 
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and 

^-aa^eloge ^ 5F = a^'^ , (2.9) 

where 

271 Jo IJ? 

Notice that we express the result (|2.71 ) directly in terms of the modification to the 
strong coupling itself, rather than that in the derived quantity cteff which was used in some 



A^TT ^ 5o.s{^^') ■ (2.10) 



previous publications [g, 0, |T 

«eff(/^^) = = - Y + • ■ ■ • (2.11) 

Although Us and ctefr are similar in the perturbative region, they differ substantially at low 
scales, and the former probably has a simpler behaviour. For example, the even moments 
of the effective coupling modification, 

r^/x^^^M/.^), (2.12) 

2tx Jo /i^ 

have to vanish for all integer values of p, whereas those of 6a g do not. The translation 
dictionary for the moments is in any case rather simple: 

^2p+i = {-irip + DnA^p+i , A2, = (-l)>A'2p , (2.13) 

where 

A, ^ ^A,, = ^ /i^^ log/.^ 5«efr(/i^) . (2.14) 

^ dv 271 Jo u 



dp 271 Jo yU^ 

Studies of the non-singlet contribution to DIS structure functions suggest that A2 



-A'2 ~ 0.2GeV2 P 



This is the formulation of the dispersive approach to power behaved corrections used 
where there is a single gluon propagator. However, the calculation which follows in section 
^ involves two gluons, and so the above argument needs to be generalised. 

Where a squared Feynman diagram contains two or more gluons, the generalisation 
of the above requires a separate independent dispersion relation and dispersive variable 
associated with each gluon. This leads to a characteristic function J?-'(ei, ■■■,€„) in all of 
these variables, and, following the above argument, we will require those terms of JF which 
are non-analytic at in all the arguments. However, in the case where there are two or 
more internal gluons all constrained to have the same 4-momentum, we can simplify this 
to require only one dispersive variable for all of them, as discussed below. 

Suppose a squared Feynman diagram contains n such gluons whose 4-momenta are all 
constrained to be equal to k, for some k. (Although only the case ri = 2 is relevant to 
phenomenology, its treatment takes no fewer lines than that of the general case.) Then, 
by defining p = —k^/Q^, we see that the dependence of the characteristic function on 
ei, ■ ■ ■ ,en is given by 

dpfip) 



(p + ei)---(p + e„) 
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where the integration hmits and the function / depend on the particular calculation. This 
may be expressed in partial fractions in the form 



En 



1 



i=l 



dpfjp) 



which yields the result 



(2.16) 



En 



i=l ^jjki 



where 



\n-l 



dpfjp) 
P + e 



(2.17) 



(2.18) 



Now the generalisation of the one-gluon result ( |2.4D to the multiple-gluon case (which 
can easily be seen from a generalisation of the argument in ^) is 



F=(-i)"n 



/i| 



p.(p|)E(-i: 



i,=0 



(2.19) 



Substituting ( ^.17] ) into this gives 



F=(-i)"n 



dfx^j 



p] 



p.(p?)E(-ip E n 



-ikpl 



i,=0 



Htkpi/Q') ■ (2.20) 



Rotating the integration contour over /i^, as in the single-gluon case, and noting that only 
terms with k = 1 and ii = 1 contribute, we obtain 



-1)' 



/if 



\n-l 



pI 



as{p\] 



as{pl] 



p] 



where 



n 

1 



i,=0 



ijp] + pi 



dp 3 I 2^ 



2m 



:DiscjF(— e) , 



(2.21) 



(2.22) 



and Joo is some integral along a contour at infinity. 
Using (|2.1|) n — 1 times one finds 



p"^ 



[asip'TGip'/Q') + /ooN . 



(2.23) 



As before, we expect to recover non-perturbative contributions to F from the modifi- 
cation in the strong coupling, 6as- The behaviour of any power correction is given by the 
non-analytic parts of J^{e), in precisely the same way as that given by ^(e) in the single- 
gluon case, but the coefficient multiplying the correction is different. In the case n = 2, we 
obtain the following non-perturbative contribution to the observable F: 



6F 



p"^ 



{2as{p')6as{p') - [5«s(p')]') G{pVQ') 



(2.24) 



We shall see that the relevant non-analytic contributions are 

^^aieloge SF = , (2.25) 

and 

^-iaselog^e ^ 5F = a^^log^ , (2.26) 
where Di and D2 are defined by: 

D, ^ r^fi' (2as{fi')6as{fi')-[6as{fi')]') , (2.27) 

logDa = r^fiHogfi' (2a,{fi')6a.{fi')-[6as{fi')]') . (2.28) 

While we expect the form of ^^(yU^), and hence Di and D2, to be universal, we have as 
yet no numerical values for them, (unlike for ^2)- It will be necessary therefore to extract 
values for Di and D2, either from experimental results or from some model of the form of 
as(/i^) (of which various models have been proposed 0). 



3 DIS Structure Functions 

We consider the deep inelastic scattering of a lepton of momentum / from a nucleon of 
momentum P, with momentum transfer q. The main kinematic variables are = —q"^, 
the Bjorken variable x = Q^/2P ■ q and y = P ■ q/ P ■ I Q"^ /xs, s being the total cm. 
energy squared. 

The differential cross section is 

_^.!p!{[i + a-,)^]f,(.) + 2(i-,)f.M} (3,1) 

where Ft{x) = 2Fi{x) and Fi{x) = F2{x)/x — 2Fi{x) are the transverse and longitudinal 
structure functions, which also have a weak dependence which we do not show explicitly. 
(For simplicity we are neglecting any contribution from weak interactions, i.e. Z" or W^ 
exchange.) 

In the parton model, to order a^, we have 



Ft{x) = (3-2) 
q 

Fl{x) = (3.3) 



where q{x) and q{x) are the quark and antiquark distributions in the target nucleon. 
The 0{as) contributions, such as those shown in figure |l], are given by 

= ^T.^l[ j{CFC.MWIi) + <l{^/0]+TRCU09{^/i)} , (3.4) 

where g{x) is the gluon distribution, Cp = 4/3, Tr = 1/2. The coefficient functions 
Cjj (0 &re the integrals in the appropriate factorisation scheme, over the final-state variable 
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7] = P- r/P • g (0 < < 1), of |T9l: 



(1-0(1-^) 



CUi.ri) = 8^(1-0- 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



m 



r]{l-r]) 

Here we concentrate on the singlet contribution, i.e. the corrections to Ci^g{^), as shown 
figure |I|. 




Figure 1: Flavour Singlet Contribution to Deep Inelastic Scattering 



4 Power Corrections in Flavour Singlet DIS 

In the normal perturbative treatment of DIS, the asymptotic freedom of QCD enables 
us to treat the initial state partons as free particles confined within the nucleon; and so 
in a singlet calculation we would start from a free gluon and convolute the perturbative 
result with the gluon distribution function g{x). We do not know how to do this in a 
calculation of power corrections, since the models we use consider modifications to the 
gluon propagator (loop insertions in the renormalon model, or, equivalently, a 'mass' in the 
dispersive approach). Let us therefore perform the calculation by considering our initial 
state gluon to be radiated from a fermionic parton. We may then try to recover the singlet 
contribution to the power corrections by deconvoluting the result with the quark to gluon 
splitting function, as performed in |]rB|, or we may leave the result as it is and interpret it 



as a second order non-singlet contribution. We might hope that these two interpretations 
would give similar predictions for power-suppressed corrections. 
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Figure 2: Diagrams Generating Flavour Singlet Contribution 



Consider the contribution to DIS from the diagrams in figure ^j. It is convenient to work 
in the Breit frame of reference which is the rest-frame of 2xP + q. In this frame 

the momentum transfer q is purely spacelike, and we choose to align it along the +z axis. 

The momentum of the initial state parton isp = xP/C,, (x < ^ < 1); and let us introduce 
the variables p = — s^/Q^, r] = P- r/P-q, f] = P- k/P-q,x the azimuthal angle between 
r and s, and 7 the azimuthal angle between r and k. There is also an overall azimuthal 
angle 0, the dependence on which is trivial. 

In the Breit frame the kinematics are given by: 

P = |g(l/x,0,0,-l/x) (4.1) 

p = ig(i/e,o, 0,-1/0 (4.2) 

q = |Q(0,0,0,2) (4.3) 

s = |Q(so,s_lcosx, s±sinx,S3) (4.4) 

r = lQ{zo,z±,0,Z3) (4.5) 

k = |Q(^o,z_LCOS7,2_Lsin7,Z3) . (4.6) 

The definitions of p, 77 and f] along with the on-shell conditions for the outgoing particles 
require that 

.„^i-P«-4 .3^-i-P« + 4 (4.7) 

2 2 

zo = V+"^ .3=r^-^ (4.8) 

-^o = V+"^ .-3=r^-|. (4.9) 

Conservation of the 0th and 3rd components of 4-momentum give the conditions 

V + V + P^ = 1 (4.10) 
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— + — + — 

Ar] 4:7] 4p^ 



(4.11) 



while conservation of transverse momentum requires that s±, z± and z±_ satisfy the triangle 
inequalities 



\z± — s±\ < z± 
\z± — z±\ < s± 



(4.12) 
(4.13) 



Two variables, a and /?, are required to parametrise the permitted values of s±, z±_ and 
z^. Let us choose to write: 



4ar/(l-r/)(l-0 

4r/(l - 

4p(l - 
1 — 



'1 — a)p^ + a?]?] — 2 cos/3A/a(l — a)ririp^ 



[1 — a)ri + ar^p^ + 2cos(3Ja{l — a)ririp^ 



(4.14) 
(4.15) 
(4.16) 



Given s±, z±_ and z±, the angles x and 7 are determined up to a sign. We may then 
choose a, P, 1] and p as the independent variables, with phase space 



< a < 1, < /? < vr, < < 1, < p < (1 - r^)/^. 



(4.17) 



To estimate power corrections to perturbative calculations, we must perform the calcu- 
lations as though the gluons had small non-zero masses pf = eiQ^ and p^ = ^2Q^- Using 
the machinery of section ^ we may write 



txTi{x\ ei) - t2Ti{x\ 62) 



ei - €2 



where 



Ti(x\ e) 



(27r)2 V ' 



a(e;e)g(a;/0 



(4.18) 



(4.19) 



The quantities Ci are the integrated perturbative matrix elements with modified gluon 
propagators. (Note that these are not the same quantities as those introduced in equation 
( p.4| ), since they refer to different diagrams.) 

To integrate the matrix elements we apply the operator 



rf3p' 



(27r)32r0 (27r)32A;0 (27r)32p'0 
Integrating away p' and making substitutions for r and k gives 

Q2 /• z±dz±d(f)dri z±dz±d'ydf] S(^A + y/B + 2z±z±cos'y 



(4.20) 



32(27r)^ 



V 



y/B + 2z±Z^_ COS 7 



(4.21) 



where A and B do not depend on 7. 
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Next wc integrate over 7. There are two values satisfying the integration condition, and 
they differ by a sign. This gives 



f dr] dr] dz^dz^ 



32(27r)^ J T] J] 2^_|_^_|_| sin7| 
Applying the parametrisation in terms of a and (3, we find that 



(4.22) 



= 32{l-0'smP^a{l-a)r,^tp/e (4.23) 



2;2^^x|sin7| = 8{l-^)smp^a{l-a)r]f]p/^ . (4.24) 
Therefore the integral operator is 



8(27r)5 

and hence 



dpdr)dad(3d4> ; (4-25) 



e 8(27r)^ Jo Jo Jo Jo Jo p + e 

where the function / is obtained from the appropriate matrix elements, excluding the 
denominators of the gluon propagators. 
The integration over is trivial. 

The part of this that is nonanalytic as e — > is the contribution to the integral near 
p = 0. Therefore we try to proceed by expressing / as a series expansion in pa about 0, up 
to C(p^). In order to produce valid series expansions, it transpires that we have to divide 
the 7] integral into three regions: 

(i) < 77 < p/k, 

(ii) p/k < T] < 1 — p^ — p/ K, and 

(iii) I- p^ - p/k <r] <1 - p^, 

where k. is some small arbitrary quantity on which the final answer should not depend. (In 
practice we perform the calculation in the limit of small k.) 

For region (ii) we may directly expand in p without any problems. The integral over (3 is 
then straightforward, being a sum of terms of the form const. cos"/9, and upon performing 
this integral the terms in non-integer powers of p vanish. The integrals over a, then 1], then 
p are not difficult, provided we note that the non-analytic parts of the p integrals are given 
by: 



^ dp (-l)"-'e"loge, (4.27) 



i 



p + e 

dp i(-l)"-^e"log'e. (4.28) 



p + e 

Note that because of the limits we have chosen for the rj integral, contributions up to C(p^) 
will also arise from all higher terms in the series expansion of /. These contributions vanish 
as K — > 0, so the calculation performed in this limit is valid. 
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For region (i) such a naive power series becomes invalid, but we may write i] = Xp and 
then expand as before. In contrast to the situation in region (ii) there are now no terms in 
non-integer powers of p. The integration over f3 may be performed by writing t = tan(/5/2), 
and that over a by making a series of substitutions. The integration over A (i.e. rj) can 
then be performed, and, since we are only concerned about the limit k ^ 0, we neglect all 
terms that vanish in this limit. The p integration then proceeds as above. 

For region (iii), we use the symmetry between the outgoing quark and antiquark: if we 
replace t] hj 1 — r] — p^ we see that the integrals over regions (i) and (iii) are equal. 

Putting all this together we have: 

^Ct = -|(2-63e + 63^2_2^3^121oge-27eioge-27^2i^g^_^12^3i^g^)lQg^ 
+1 (4 + 3e - 3^' - 4^=^ + log e + log (log e - 1 + I log e) log e 



+ 1(2 + 25^' - 25^' - 2^^ + ISe'loge + 15^=^ log doge 



-2(5^' - 5^' + 2e\og^ + 2edogO(loge - 1 + iloge)eloge , (4.29) 

3* 



^Cl = -|(l-3e + 2e=^-3e'log01oge 



-215(17 + 75^2 - 125^^ + 33^^ + 301og^ + 75^^ log^ - 45^^1ogOeloge 
+^(2 - 15^2 + lOe + 3^^ - 15e'logO(loge - 1 + iloge)eloge . (4.30) 

The terms that diverge as e — are responsible for the logarithmic scaling violations to 
the structure function (c.f. [^). To see this, note that the leading divergent pieces (i.e. those 
proportional to | log^ e for Ct and to loge for Cl) can be expressed as a convolution of the 
quark to gluon splitting function with the first order gluon coefficient functions: 

2(4 + 3e-3^--4f^+6^1ogt + 6eMosO ^ 2r + (1 - C)> ^ + (4,31) 

^'^^^^ + ^1'-^^'^'°^^' ^ 8J(l-0«i±il^. (4.32) 

The subleading divergence in Ct, i.e. that proportional to loge, is factorisation-scheme- 
dependent. 

The terms proportional to e log e and e log^ e give the power corrections. These 

are consistent with the results in []T^. Comparison of the results presented here with those 
T8|1 indicate that the terms containing the factor (log ^ — 1 + | log e) correspond to 



m 



the contributions from the double renormalon pole, while those terms without this factor 
correspond to the contributions from the single renormalon pole. These results also suggest 
signs for the different contributions: the contributions in |jl8|] from the single poles appear 
here with a positive sign whereas those from the double poles appear with a negative sign. 
This leads to power corrections to the structure functions given by: 

= ^ E 4 f T^CM)q{x/0 , (4.33) 



where 



■2(2 + 25^2 _ 25^3 _ 2^5 ^ ^5^2 ^ + 15^3 j^g ^) 
. 5^ 
-2(5^ - 5e + 2^ loge + 2e logO log ^ 



(4.34) 
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8(17 + 75^2 _ 125^3 ^ 33^5 ^ 33 i^g^ + 75^3 j^g^ _ 45^5 j^g^) 



225^ 

, 8(2-15e^ + 10e^ + 3e^-15e=^logO , 
+ — log 



(4.35) 



5 Results and Conclusions 



The l/Q"^ power corrections arising from the diagrams shown in figure ^ are given above in 
equations ( ^4.34 ) and ( 4.35| ). While we do not know the values of Di and D2, we can still 
make some qualitative predictions. Figure]^ shows plots of Kt{x) and Kl{x), defined by 

Di TrCf 



SFdx) 



(27r): 



-KAx) 



(5.1) 



i.e. Ki{x) are the coefficients of the l/Q"^ power correction, excluding the unknown factor 
Di. These were calculated at = 500GeV^, using the corresponding MRST (central 
gluon) parton distributions [^. The value of D2/e was set to be O.OGGeV^, i.e. approx- 
imately A^, following the result of 



ISj . (The qualitative behaviour of the Ki does not 



change provided we keep D2 <^ Q ) 




Figure 3: Graph showing Kt{x) and Kl{x). 



As can be seen from the graph, Kt^x) and Kl{x) both tend to zero at large x, and 
both diverge at small x. Kt{x) is positive and Kl{x) is negative, with the magnitude of 
Kl considerably larger at small x than that of Kt- The corresponding quantity related to 
F2/X, which is K2/X = Kt + Kl, also therefore behaves qualitatively like Kl. 

Therefore if Di is positive, we have a negative correction to both F2/X and Fl. If Di 
is negative then the correction is positive. In either case the corrections are important at 
small X. 

We do not know the sign or magnitude of Di, and there is no a priori reason why it 
should be either positive or negative. However, since A2 is positive, we presumably have 



11 



Sasin"^) being predominantly positive. We therefore might conjecture that Di would also 
be positive. 
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